A new determination of the electromagnetic nucleon form factors from QCD Sum 

Rules 
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We obtain the electromagnetic form factors of the nucleon, in the space-like region, using three- 
point function Finite Energy QCD Sum Rules. The QCD calculation is performed to leading order 
in perturbation theory in the chiral limit, and also to leading order in the non-perturbative power 
corrections. For the Dirac form factor, Fi(q 2 ), we get a very good agreement with the data for 
both the proton and the neutron, in the currently accessible experimental region of momentum 
transfers. Unfortunately this is not the case, though, for the Pauli form factor F2(q 2 ), which has a 
soft q 2 -dependence proportional to the quark condensate < 0|qq|0 >. 



The determination of the electromagnetic nucleon form 
factors is an old standing problem in QCD. For a re- 
view, see pj. Calculations based on perturbative QCD 
(PQCD), together with sum rules estimates for the nu- 
cleon wave function, are difficult to compare with data 
due to the extreme asymptotic nature of these theoreti- 
cal results. Recently, a new analysis based on light-cone 
QCD Sum Rules [J has been carried out improving the 
agreement with data from within a factor 5-6 to a fac- 
tor of two. Here we attempt a Finite Energy QCD Sum 
Rules (FESR) determination of the Dirac Fi(Q 2 ) and 
of the Pauli ^(Q 2 ) form factors, in the region of ex- 
perimentally accessible momentum transfers. The QCD- 
FESR approach is interesting, since power corrections as- 
sociated to vacuum condensates of different dimensions 
decouple at leading order in PQCD. 

As it is well known this technique is based on the Oper- 
ator Product Expansion (OPE) of current correlators and 
on the notion of quark- hadron duality . Our calcula- 
tion will be done to leading order in PQCD, in the chiral 
limit including also the leading-order non-perturbative 
power corrections associated to the quark-condensate and 
to the four-quark condensate. 

By considering the interpolating current with proton 
quantum numbers 

Vn(x) = e abc [u a {x){C la )u b {x)] ( 7 Vd c (x)), (1) 

and the electromagnetic current 
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we are interested in the three-point correlator 
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-ip' - p) 2 > is fixed. See Fig.l. 



where Q 2 = —q 2 

The current Eq.l couples to a nucleon of momentum p 
and polarization s according to 



where u(p,s) is a nucleon spinor and Ajv is a phenomeno- 
logical parameter that gives us the current-nucleon cou- 
pling. This parameter has been estimated, for example, 
using two-point QCD sum rules involving the current rj^ 
i-fl. 

Going to the hadronic sector, after inserting a one- 
particle nucleon state, the three-point function (3) can 
be written in terms of the nucleon form factors Fi(Q 2 ) 
and F2(Q 2 ), defined as 



(fcisi |J^ M (0)|fc 2 s 2 ) =«jv(*i,*i) 
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u N (k 2 ,s 2 ) , (5) 



where q 2 = (k 2 — fci) 2 , and k is the anomalous magnetic 
moment in units of nuclear magnetons (n p = 1.79 for 
the proton, and n n — —1.91 for the neutron). The form 
factors Fi^{q 2 ) are related to the electric and magnetic 
(Sachs) form factors G_e(q 2 ), and Gm(q 2 ), measured in 
elastic electron-proton scattering experiments, according 
to 



G E (q 2 ) = F^q 2 ) 
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G M (q 2 ) = F^q 2 ) + KF 2 {q 2 ) , 



F2(q 2 ) , (6) 
(7) 



(0 1^(0)1 N(p,s))=\ N u(p,s) 



where G P E (0) — 1, G P AI (0) = 1 + k p for the proton, and 
G^(0) = 0, G^(0) = K n for the neutron. Next we com- 
pute the hadronic spectral function by inserting a com- 
plete set of intermediate nucleonic states in (3) and com- 
(3) puting the double discontinuity in the complex p 2 = s, 
p' 2 = s' plane. If we stay with s, s' < 2.1GeV 2 , i.e. below 
the Roper resonance, we can approximate the hadronic 
spectral function by the single-particle nucleon pole plus 
a continuum with thresholds sq and s' (sq,s' q > M^) 
that we expect will coincide with the PQCD spectral 
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FIG. 1: The three-point function, eq. 3, to leading order in 
PQCD 



function (local duality). In this way we get 

Imn M (s, s', Q 2 )\had = 7T 2 X 2 N 6{s - M 2 N ) 
xd(s> - MDiF^A, + ^ [ -F 2 (q 2 )B, v 



q v }@ (s - s) 



Imll^s, s',Q 2 ) 



PQCD 



e(s-so), (8) 



where for simplicity we set sq = s and and B^ v 
correspond to the following tensor structures 



and 



B^ v = i>cF^ v i> + M N (fia,j, v + o^f) 
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Going to the QCD sector, to leading order in PQCD 
and in the chiral limit, we have to calculate the imagi- 
nary part of the diagram shown in Fig.l. The important 
point is that there are several Lorentz structures, anal- 
ogous to those we found in the hadronic sector. Before 
invoking local duality it is necessary to choose a partic- 
ular Lorentz structure present both in the QCD as well 
as in the hadronic sectors. 

The term f'^i^'p turns out to be appropriate. It allows 
to project Fi(Q 2 ) since this structure does not appear 
together with F 2 (Q 2 ) in the hadronic spectral function. 
On the other hand, due to vanishing traces, the quark 
condensate to be considered later, also does not involve 
this structure. In principle, however, there are four quark 
condensate terms associated with such structure. How- 
ever, those terms do not contribute to the FESR since 
the associate double discontinuity vanishes. After a very 
lengthy calculation, the imaginary part of the perturba- 
tive expression of the correlator, associated to the desired 
structure /5'7 M ^, can be written as 



ImIP(s, s', Q 2 ) = frfi (a + /3[-323 Q 12 - 
Q w Pio(s, s') - 10 Q 8 P 8 (s, s') + Q 6 P 6 (s, s') 
+ Q A Pi{s, s') - Q 2 P 2 (s 1 s') - P {s, s')]) + ... (11) 

In the equation above the dots denote the terms associ- 
ated to other Lorentz structures and we have introduced 

/j- 1 =46087r 2 [ g 4 +(s - s') 2 + 2Q 2 (s + s')] 1 
323Q 2 + 378 (s - s') 
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4608tt 2 , 



(12) 



and the set of polynomials Pi(s, s') given by 

Pio = 1993s + 1237s', 

P s = 512s 2 + 323ss' + 134s' 2 , 

P 6 = -7010s 3 + 1188ss' 2 + 550s' 3 , 

Pi = -5395s 4 + 7010s 3 s' + 2610s 2 s' 2 + 3146ss' 3 + 2165s' 4 

P 2 = (s - s') 2 (2213s 3 - 2589s 2 s' - 3099ss' 2 - 1567s' 3 ) 

P = -378(s - s') 4 (s 2 - 2ss' - s' 2 ) 

It is interesting to mention that we got both explicit 
terms, where the desired tensor structure was there from 
the beginning, as well as implicit terms, i.e. those terms 
where the tensor structure emerged only after performing 
the integrals. The next step is to invoke global quark- 
hadron duality in the frame of the FESR, which, as we 
mentioned, are organized according to dimensionality. 
The FESR of leading dimensionality are 

so pSq — s 

ds / ds'Imn(s,s',Q 2 ) \ HAD = 
ii Jo 

so rs —s 

ds / ds'ImII(s,s',Q 2 ) \ QCD ■ (13) 
o Jo 

We have chosen a triangular region to integrate in the 
s, s' plane, but the result is quite independent from the 
integration region [|| and 0- In this way one obtains 



FtiQ 2 ) 



9216 vr 4 (Q 2 + 2s )Ajv 2 



x (A + B ln( 
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Q 2 + 2s 
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(14) 



where we have defined 

A = 2 s (96 Q 6 + 297 Q 4 s + 158 Q 2 s Q 2 - 112 s 3 ) and 

B = 3(Q 2 + 2 s ) (32 Q 6 + 67Q 4 s + 7Q 2 s 2 - 21s 3 ). 

Notice that in the previous equation we have the stan- 
dard logarithmic singularity arising from the chiral limit. 
The leading asymptotic term turns out to be 



lim Q 4 F^Q 2 
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2560 7T 4 X 2 N 



(15) 



Qualitatively, this asymptotic behaviour agrees with 
expectations. From QCD sum rules for two-point func- 
tions involving the nucleon current (1) it has been found 



The problem is that the asymptotic behavior does not 



G En (Q 2 ) 



An = 0.01 * GeV^ 
s = 1-2GeV2 



Q2 (GeV^) 



FIG. 2: Theoretical results (solid line) versus corrected ex- 
perimental [3 data on Fi(Q 2 ) 
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(1 - 3) x l(T 2 GeV 3 , and ^ 



(1.1 — 1.5) GeV. The higher values of Xn and so come 
from Laplace sum rules Q, and the lower values are 
from a FESR analysis p| which yields the relation s|] = 
192tt 4 A^. After fitting Eq.(14) to the experimental 
data, as corrected in 0, we find A A? = O.OllGeV 3 , and 
So = 1.2GeV 2 , in line with the values discussed above. 
Numerically, so is well below the Roper resonance peak, 
thus justifying the model used for the hadronic spectral 
function. The predicted form factor F\(q 2 ) is shown in 
Fig. 2 (solid line) together with the data. The agreement 
is quite impressive. 

From the two leading power corrections in the OPE, 
without gluon exchange, the one proportional to the 
quark condensate does not contribute to Fi(Q 2 ), while 
the other, proportional to the four-quark condensate, has 
a vanishing double discontinuity in the (s, s') complex 
plane. For details see the original article [lCj . 

In order to extract F 2 (Q 2 ), we have to consider the 
leading-order non-perturbative corrections to the OPE, 
which in this case corresponds to the quark condensate. 

In the case of the proton, the contribution involving the 
up-quark condensate vanishes (due to vanishing traces) 
and therefore we only have a piece proportional to (dd) . 
Our choice of Lorentz structure in this case is ^7 M , which 
appears in the QCD sector as well as in the hadronic 
sector multiplying FAQ 2 ) but not Fi(Q 2 ). We refer to 
the original article [10( for the full expressions. Here we 
will give only the final result that emerges form the FESR 
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FIG. 3: Experimental data on Gb(Q 2 ) for the neutron [Tll |. 
together with the theoretical results 



agree with the expectations. We find 



lim F 2 (Q 2 ) 

Q->—>oo 



(dd) 



18k p M n tt 2 X 2 n 



(17) 



and we would expect F 2 (Q 2 ) to fall faster than Fi(Q 2 ) 
at least by a factor of X/Q ||- A comparison of F 2 (Q 2 ) 
from equation (16) with data shows a disagreement at the 
level of a factor two, which cannot be improved adjusting 
the values of A and So- The main reason behind the 
disagreement is the soft <? 2 -dependence of F 2 (Q 2 ). 

We can do the same analysis for the neutron form fac- 
tors. It turns out that Fi n (Q 2 ) for the neutron is numer- 
ically very small and consistent with zero, except near 
Q 2 = due to the divergence in the chiral limit. Since 
Fin(Q 2 ) ~ 0, the Sachs form factor is proportional to 
F 2n (Q 2 ). In Fig. 3 we show a plot of the electric Sachs 
form factor for the neutron. At low Q 2 there is reason- 
able agreement with the experimental data. However, for 
higher momentum transfers the disagreement turns out 
to be serious due to the soft l/Q 2 behavior of F 2n (Q 2 ). 
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